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ABSTRACT 

O ! Using the results of recent work in shear instabilities in stratified fiuids, we 

^ ! show that the resonant interaction between large-scale fiows in the accreted H/He 

^ . envelope of white dwarf stars and interfacial gravity waves can mix with the star's 

Q ! envelope with the white dwarf's surface material, leading to the enhancement 

of the envelope's C/O abundance to levels required by extant models for nova 
i outbursts. 



Subject headings: stars: novae 



> 
OO 

O I The substantial enrichment of CNO nuclei in the ejecta from novae (cf. Truran 1985; 

Gehrz et al. 1998 and references therein) has been a puzzle for over two decades. Early 
theoretical models of nova outbursts (e.g., Starrfield, Truran & Sparks 1978; Fujimoto 1982) 
^ ■ clearly showed that nuclear processing during hydrogen burning during the nova fiash could 
I ' not account for the observed CNO abundances, which can reach 30% by mass. These early 
' studies already recognized that the solution to the puzzle must involve "dredgeup" of C/O 



^ ■ from the white dwarf before or during the nova outburst. This mixing was required both to 
meet the constraints on CNO abundances in the ejecta and to power the nova itself, since 
the energy production rate per unit mass depends directly on the metallicity (e.g., Wallace 



' & Woosley 1981); thus, Starrfield et al. (1978) and Fujimoto (1982) showed explicitly with 
1-D models that runaway in a pure H/He envelope did not release enough energy in order 
to eject enough matter with sufficient velocity to match observations. 

Concerns that mixing may occur at the interface between the accreting matter and the 
underlying star (and already accreted material) followed closely upon the recognition that 
such mixing was essential in order to understand the elemental composition of the nova 
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ejecta (Starrfield et al. 1972). At that time, there was aheady some interest in understand- 
ing mixing at the interface between a stellar surface and an accreting flow. For example, 
Kippenhahn & Thomas (1978) examined shear flow instability in the stratified boundary 
layer between a white dwarf and the infalling accretion flow associated with an accretion 
disk; and established the linear stability properties (based on using the Richardson number^ 
as the control parameter). 

The shear instability considered by Kippenhahn & Thomas (also Sung 1974) has been 
extensively revisited (viz., MacDonald 1983). Kippenhahn & Thomas conjectured that this 
instability saturates at the marginal state for stability, and therefore weak mixing; Mac- 
Donald, upon revisiting this problem, showed that the shear instability would lead to rapid 
dispersal of the accreted matter over the entire white dwarf surface (as opposed to the rela- 
tively narrow accretion belt which emerged from Kippenhahn & Thomas' analysis), but also 
suggested that the radial mixing time was long (set by the thermal timescale of the enve- 
lope). These arguments lead to rather minimal mixing, and for these reasons shear mixing 
has not been regarded as a likely candidate for the required mixing process.^ Indeed, until 
the late 1990's the absence of a plausible mixing mechanism was considered to be a major 
stumbling block for understanding novae. In the mid-90's, several authors conjectured that 
the convection which was known to initiate some ~ 1,000 years before runaway might be as- 
sociated with convective undershoot and convective penetration, processes which might lead 
to mixing of the stellar C/0 into the envelope (Shankar, Arnett, & Pryxell 1992; Shankar 
& Arnett 1994), but quantitative calculations were not done until the mid and late 1990's 
(Glasner & Livne 1995; Glasner et al. 1997; Kercek, Hillebrandt, & Truran 1998a). These 
more recent calculations investigated the possibility that convective undershoot just before, 
and possibly during, nova runaway might lead to the required mixing. However, Kercek, 
Hillebrandt, & Truran (1998b, 1999) have shown convincingly (both by comparing two and 
three-dimensional simulations, and by conducting resolution studies in which the extent of 
mixing was measured as a function of grid resolution) that convective undershoot was not 
likely to work as an effective mixing process. In particular, the resolution studies showed less 
mixing as grid resolution was increased. This can be readily understood if the boundary layer 



^The Richardson number, Ri = agAlo/U^, is a measure of the competition between the stabiUzing effect 
of buoyancy and the destabihzing effect of the shear flow. Here a is the coefficient of volume expansion, 
g is the local gravitational acceleration, lo is a local characteristic length scale, and Uo is the shear flow 
amplitude; A is the Atwood number across the material interface separating the white dwarf surface and the 

H/He envelope [= (pstar - Penvelope)/(Pstar + Penvelope)]- 

^However, very recently, Briiggen & Hillebrandt (2001a,b) have begun to re-examine the nonlinear as- 
pects of this problem computationally, in an attempt to place the earlier analytical calculations on a firmer 
quantitative footing. 
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between the stellar surface and the accreted (convecting) envelope is laminar: in that case, 
since the dominant viscosity in these simulations is numerical, increased resolution leads to 
a thinner boundary (or mixing) layer, and whence to less mixing as the grid resolution is 
increased. Thus, it would appear that we are once again lacking an effective mixing process. 

For this reason, we have recently reexamined the physics of shear flow instabilities 
(Alexakis, Young, & Rosner 2001; Young ct al. 2001). The question to answer was whether 
previous astrophysical studies of this subject had in fact fully explored this mixing process. 
As we show below, the past work in fact missed an important aspect of shear mixing in 
stratified media. In the following, wc will extract the critical aspects of our earlier results 
that apply to the problem at hand. 

One very important aspect of the shear mixing problem is that the instability that 
leads to turbulent mixing between fluid layers depends only on certain essential features of 
the shear flow. For present purposes, it suffices to consider prototypical velocity profiles 
of the form U{z) = Uo + Ui\n{z/a + 1) for z > {U{z) — below the interface), or 
U (z) —Uo + Ui tanh.{z/a) for 2; > {U{z) — below the interface), where [/„ is the velocity 
jump (if any) at the envelope/stellar surface interface, z is the vertical coordinate (with 
z — marking the initial envelope-star interface), and a is the characteristic scale length of 
the shear flow in the envelope.^ It is then well-known that if Uo = 0, the Kclvin-Helmholtz 
instability is entirely absent for a velocity profile with either of the two functional forms 
given above. Since (in the presence of viscosity) the relative velocity between the two layers 
of fluid at the interface must be zero (i.e., an attached flow boundary condition) and since 
a is proportional to the viscous boundary layer thickness, Kelvin-Helmholtz instability is 
unlikely to be important in contributing to the mixing on relevant spatial and time scales 
for the above wind profile, which appears to be a reasonable approximation to the actual 
boundary layer flow (cf., Alexakis et al. 2001). 

Now, the same argument applies to the generation of terrestrial surface water waves by 
winds; and as originally pointed out by Miles (1957; see also Phillips 1957), winds are nev- 
ertheless able to amplify such surface waves to finite amplitude. Hence, there must be some 
other instability present. One such instability (critical-layer instability) was identified and 
studied extensively in the linear regime by Miles (1957), Howard (1961), Lighthill (1962), and 
others. This instability originates from the continuum of unstable modes formed when sur- 
face gravity waves travel at the same velocity as the wind at some height above the interface. 



^The logarithmic velocity profile is commonly observed in the boundary layer of winds blowing over the 
surface of extensive bodies of water [cf. Miles 1957]; the tanh profile has the advantage of bounded shear 
velocity far from the shear interface. 
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In order to treat this instability, it is essential to formulate the shear instability problem more 
generally; this is in part the motivation for considering velocity profiles of the forms given 
above; velocity profiles of this type in stratified atmospheres have been explored extensively 
by the geophysical fluid dynamics community. These previous geophysically-motivated stud- 
ies were largely confined to the parameter regime characteristic of the water/air interface; 
but recently Alexakis, Young, & Rosner (2001) have fully explored the control parameter 
space governing these instabilities: these nondimensional parameters are the Atwood num- 
ber A and the Richardson number Ri (see Alexakis et al. 2001); in physical terms, the key 
parameters are the gravitational acceleration g, the shear scale length a, the density ratio 
Pstar/Penveiope; ^-^d the shear amplitude U. In this Letter, we now apply the results presented 
by Alexakis et al. to the astrophysical context by developing a new model for the interface 
mixing. 

First, consider the underlying physics: This is most readily done in the context of a 
particularly simple model for the interface, in which the shear flow has step discontinuity 
across the density interface between the C/0 white dwarf surface and the bottom of the H/He 
envelope. Start with the simplest case, in which we ignore stratification. The classical Kelvin- 
Helmholtz instability is then based on the observation that a spatial interface perturbation 
can be destabilized because the flow must speed up over the "hills" of the perturbation and 
slow over the "valleys" ; BernouUi's law then tells us that a low-pressure region develops over 
the "hills", and a high-pressure region over the "valleys", thus puUing up the "hills" and 
pushing down the "valleys" , leading to a linear instability whose growth rate 7 ~ kU, where 
k is the wavenumber of the interface perturbation and U is the shear amplitude. 

If the shear flow interface is not a step, but has flnite thickness (viz., given by a, as 
above), and if stratiflcation is allowed, then it is well-known that the dispersion relation 
is no longer linear (Chandrasekhar 1962, §102), and that both low and high wavenumber 
cutoffs appear, with 7 > only for (^min < ka < Cmax (where the values of C^^ain and Cmax 
depend on the speciflcs of the velocity proflle (cf. Figures 119 and 120 in Chandrasekhar 
1962), and maximum growth at, for example, {ka)"^ ~ 0.5 for the tanh(2;/o') shear profile. 
As a result, instability can only occur in a finite region of the wavenumber- Richardson 
number plane; for the tanh velocity profile, the stabihty boundary is defined by the curves 
J = and J = {kaY{l — (ka)"^), with instability only in the domain bounded by these two 
curves. In order to apply this to the nova case, we simplify the actual case by assuming 
an exponentially-stratified background atmosphere of the form p{z) ~ poexp(— /Sz), with a 
horizontal shear layer of form U — UotBnh{z/a) located at the white dwarf surface; it is 
readily seen that in the event that this surface shear fiow is driven by thermal convection 
in the overlying envelope, then the unstable modes will lie in a wavelength band defined by 
6 X 10^ cm < Aunstabie < 2 X 10^ (T/10* K)^/^ {Uo/10^ cm s"^) cm, where we have assumed 
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a white dwarf of radius ~ 10~^i?sun, a shear layer thickness a ~ 10^ cm, gravitational 
acceleration g^^ ~ 2.7 x 10^ cm s~^, and density scale height /?~^ ~ 3 x 10^(T/10^ K) 
cm; inclusion of the density jump at the (C,0/H,He) interface would lower the wavelength 
of unstable modes yet further. The upper bound on this mixing scale is of order the grid 
resolution in the currently highest-resolution calculations (viz., Glasner, Livne & Truran 
1997; Kercek, Hillebrandt & Truran 1999), consistent with the observation by these authors 
that little shear mixing occurred in their computations. Since there is an upper bound on 
the shear flow length scale a in order for Kelvin- Helmholtz instability to occur at alF, and 
since the mixing scale is at most of order 10 times the shear scale, this suggests that Kelvin- 
Helmholtz instability will not be an effective CNO mixing process under any circumstances. 

In contrast, consider the interaction of the same wind with the normal modes supported 
by the interface between the stellar surface and the accreted envelope. These normal modes 
are akin to "deep water waves" seen at the surface of terrestrial oceans; and are known to 
grow in amplitude as a result of the resonant interaction between these waves and the wind. 
More specifically, at any given wavenumber A;, linear theory provides the wave's phase velocity 
I'phase = '^1^ ~ (^/^)^^^(fl' + T.k'^ / p"""^"^) ~ {Ag/ky/"^, where S is the surface tension; in the 
case of the gaseous media characterizing stars, the surface tension term is of course absent. 
For any given wind profile, U{z), where z is the vertical coordinate, one can then satisfy a 
resonance between the wind and a surface mode such that U{z) — Vphase ~ {Ag/ky/"^] that 
is, a wave with wavenumber k ~ Ag/U{zY will be driven resonantly unstable. (For typical 
values of a, and U characteristic of a white dwarf surface, one finds that the wavelength 
of unstable modes lies in the range 0.01 - 1 km.) The key issue is then how to determine 
the mixing layer width once these unstable modes cease their growth and finally saturate: 
naively, one might expect the saturation process to simply limit the mode amplitude, and 
thereby determine the width of the mixing layer. In the case of intcrfacial gravity modes, 
however, saturation is well-known to occur via wavcbrcaking (sec Chen, Kharif, Zalcski, 
& Li 1999 and references therein); it is the resulting spray that then determines (from a 
statistical point of view) the effective mean width of the mixing layer - this width can be 
substantially larger than the mode amplitude at saturation, as is well-known in the case of 
wind-driven spray from breaking ocean waves. In any case, let us assume for the moment that 
we have determined this layer width, which we shall denote as A. Finally, we note that while 
one would need in general to take into account stratification effects (viz., molecular weight 
gradients) on either side of the density jump, such effects are to lowest order unimportant 
here because the mixing layer defined by wave breaking is likely to be much narrower than 



^This bound is computed from the stability criterion J < 1/4 for the tanh velocity profile; thus a < 
1.67 X 10^ (T/10* Kf/'^iUo/lQ^ cm s'^) cm in order for Kelvin-Helmholtz instability to occur at all. 
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the local gravitational scale height. 

We are now ready to describe our simplified model: consider first the amount of carbon 
and oxygen in the breaking wave mixing layer, which we write in the form M^^q^ '^■^^'^ ~ 
apo'^^^^^^^AA^, where a is the coefficient for the C+0 mass fraction in the mixing layer, A 
is the mixing layer width (both a and A are to be determined from simulations; cf. Young 
et al. 2001), p^'^^^°P'^ is the density of the envelope at its base (i.e, in the breaking wave 
mixing layer), A is the characteristic length scale of the large-scale circulation (which can 
be identified with the outer scale of convection in the envelope); and ^ is a length scale 
transverse to the wind direction (this dimension will drop out of our formulation). Note 
that the remaining parameters appearing in this relation can be obtained from extant (1-D) 
nova models. Now, as mentioned earlier, the amount of C+O needed to be mixed into the 
envelope is roughly 1/3 by mass of the ejecta mass ilClf°''', or M^7g°P" ~ (l/3)7lC;f°''' ~ 
(1/3)- ^(2/3)ApQ°^^'°''^ • A^^ . The "sweepout time" , i.e., the time scale on which the boundary 
mixing layer is swept out by a penetrating convective roll, is just Tg^eep ~ A/C/, so that the 
time needed to mix the necessary amount of carbon and oxygen into the envelope is just 

M envelope // n ^-mixing layer ; \ 
C+O /(-'^^C+O /Tsweep), or 

Tmixing ~ (2/9)AVaAC/ 

~ 5a-^(A/10^cm)^(A/102cm)"^ ■ 
• (t//10^cm s~^)"^ years, 

with a ~ 0.3 — 1. Thus, it is evident that the evolution time scale for the envelope prior to 
nova nmaway (which is roughly of the order of the time between onset of envelope convection 
and runaway, or ~ 10'^ years) is much longer than the mixing timescale. This confirms that 
resonantly-driven mixing at the star-envelope boundary can be an efficient mixing process 
during the pre-nova star evolution; the clear next step is to verify these results via simulations 
of weakly compressible fiuids subject to these mixing instabilities. We also note that this 
mixing time scale is much longer than the dynamical time characteristic of the nova runaway 
itself. For this reason, the amount of additional C-l-0 material mixed in during the outburst 
itself can be regarded as a small perturbation. One remaining significant issue relates to 
the possible effects of magnetic fields on the C-l-0 mixing process; that is, one might be 
concerned that turbulent mixing may be suppressed if local magnetic fields in the envelope 
become large as convection sets on ~ 1,000 years before runaway. We arc not currently in a 
position to resolve this possible problem, but only note that because the conservative mixing 
time scale Tmixing << 1, 000 yrs, substantial mixing suppression by magnetic fields could be 
accommodated within this model without vitiating the main point, namely that resonant 
instability of the C-|-0/envelope boundary can lead to effective mixing across that boundary. 
This is a critical point for any nova model because novae have been observed for white dwarfs 
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with relatively strong magnetic fields (e.g., V1500 Cygni 1975; Stockman, Schmidt & Lamb 
1988). However, in the absence of a detailed calculation, this point remains to be addressed 
by theory. 

To conclude, by using the results of hnear stabihty theory, as well as extrapolating 
from existing numerical simulations of nova outbursts, we have estimated the mixing zone 
parameters, and have shown that pre-nova erosion of the wave-breaking mixing layer by slow 
convection could mix sufficient C/0 into the accreted H/He envelope to satisfy observations. 
Wc have constructed a simple mixing length subgrid prescription to describe this mixing 
process, and have shown that this subgrid model only needs to be used for the pre-nova phase. 
Further mixing during the outburst is no longer required. Because the C/0 abundance in 
the envelope builds up gradually during the pre-nova slow convective phase, we expect the 
nova envelope mass attained before outburst may be substantially larger than in standard 
models assuming a "pre-seeded" envelope. 

This work has been supported by the DOE- funded ASCI/Alliances Center for Astro- 
physical Thermonuclear Flashes at the University of Chicago. 

REFERENCES 

Alexakis, A., Young, Y.-N., & Rosner, R. 2001, Phys. Rev. E, submitted. 

Briiggen, M., & HiUebrandt, W. 2001a, MNRAS, 320, 73. 

Briiggen, M., & HiUebrandt, W. 2001b, MNRAS, 323, 56. 

Chandrasekhar, S. 1962, Hydrodynamics (New York: Dover). 

Chen, C, Kharif, C, Zaleski, S., & Li, J. 1999, Phys. Fluids, 11, 121. 

Fujimoto, M.Y. 1982, ApJ, 257, 752. 

Gehrz, R.D., Truran, J.W., Wilhams, R.E., & Starrfield, S. 1998, PASP, 110, 3. 

Glasner, S.A., & Livne, E. 1995, ApJ, 445, L149. 

Glasner, S.A., Livne, E., Truran, J.W. 1997, ApJ, 475, 754. 

Howard, L.N. 1961, JFM, 10, 509. 

Kercek, A., HiUebrandt, W., & Truran, J.W. 1998a, A&A, 337, 379. 

Kercek, A., HiUebrandt, W., & Truran, J.W. 1998b, in Proc. Workshop on Fluid Dynamics, 
Bordeaux, France, p. 14. 

Kercek, A., HiUebrandt, W., & Truran, J.W. 1999, A&A, 345, 831. 



-8- 



Kippenhahn, R., & Thomas, H.-C. 1978, A&A, 63, 265. 

Lighthill, M.J. 1962, JFM, 14, 385. 

MacDonald, J. 1983, ApJ, 273, 289. 

Miles, J. 1957, JFM, 3, 185. 

Miles, J. 1959a, JFM, 6, 568. 

Miles, J. 1959b, JFM, 6, 583. 

Miles, J. 1962, JFM, 13, 433. 

Miles, J. 1967, JFM, 30, 163. 

Miles, J. 1993, JFM, 256, 427. 

Miles, J., & lerley, G. 1998, JFM, 357, 21. 

Phillips, O.M. 1957, JFM, 2, 417. 

Shankar, A., & Arnett, W.D. 1994, ApJ, 433, 216. 

Shankar, A., Arnett, W.D., & Fryxell, B.A. 1992, ApJ, 394, L13. 

Starrfield, S., Sparks, W.M., & Truran, J.W. 1974, ApJS, 28, 247. 

Starrfield, S., Truran, J., & Sparks, W.M. 1978, ApJ, 226, 186. 

Starrfield, S., Truran, J., Sparks, W.M., & Kuttcr, G.G. 1972, ApJ, 176, 169. 

Stockman, H.S., Schmidt, G.D., & Lamb, D.Q. 1988, ApJ, 332, 282. 

Sung, C.-H. 1974, AA, 33, 99. 

Truran, J.W. 1985, in Nucleosynthesis, eds. W.D. Arnett & J.W. Truran (Chicago: Univ. 
of Chicago Press), p. 292. 

Wallace, R.K., & Woosley, S.E. 1981, ApJS, 45, 389. 

Woosley, S.E. 1986, in Nucleosynthesis and Chemical Evolution, ed. B. Hauck, A. Maeder, 
& G. Magnet (Sauverny: Geneva Observatory). 

Young, Y.-N., Alexakis, A., Rosner, R., et al. 2001, in preparation. 



This preprint was prepared with the AAS IM^jX macros v5.0. 



